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1  Executive  Summary 


The  central  theme  of  the  work  on  this  contract  has  been:  1)  the  development  of  techniques 
and  ideas  for  the  control  of  the  waveforms  of  single  photons  and,  2)  the  extension  of  the 
concepts  of  nonlinear  optics  into  the  x-ray  portion  of  the  electromagnetic  spectrum.  Ad¬ 
vances  in  the  first  portion  of  this  work  include:  the  first  demonstration  of  use  of  spread 
spectrum  technology  at  the  single  photon  level  (Belthangady  et  ah);  and  in  particular  the 
demonstration  that  it  is  possible  to  hide  a  single  photon  in  an  environment  of  a  much  higher 
density  of  photons  with  the  same  spectral  power  density,  and  2)  The  first  demonstration 
of  the  technique  of  chirp  and  compress,  again  at  the  level  of  a  single  photon  (Sensarn  et 
ah).  Key  contributions  and  advances  in  the  hard  x-ray  portion  of  this  work  include:  1) 
The  invention  of  a  new  technique  for  the  generation  of  polarization  entangled  x-ray  bipho¬ 
tons  (Shwartz  et  ah).  2)  The  demonstration  of  parametric  down  conversion  in  the  Langevin 
regime.  3)  The  first  demonstration  of  x-ray  optical  mixing.  (This  work  was  part  of  a  team 
led  by  Ernie  Glover.  4)  And  very  recently,  the  first  demonstration  of  second  harmonic  gener¬ 
ation  at  x-ray  wavelengths.  Additional  contributions  include  the  demonstration  of  a  greatly 
simplified  technique  for  obtaining  temporally  long  biphotons  (Chih-Sung  Chun),  as  well  as 
the  suggestion  of  a  second,  backward  wave  technique  that  has  not  yet  been  demonstrated. 

In  the  following  sections  we  follow  the  AFOSR  final  report  guideline  and  will  not  repeat 
already  published  and  well  documented  work.  For  completeness,  in  these  cases,  we  will 
include  the  published  abstract.  When  the  work  is  still  underway,  and  not  published  or  not 
completely  published,  we  expand  and  give  details  here. 

We  note  that  this  work  has  been  jointly  supported  by  the  US  Army  Research  Office. 
Before  this  contract  began,  we  were  jointly  supported  by  DARPA,  and  the  work  described 
here  benefited  from  that  interaction. 

We  begin  by  summarizing  Publications  and  Presentations  during  this  contract  period. 
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2  Publications  and  Presentations 


2.1  Publications 

(1)  Chinmay  Belthangady,  Chih-Sung  Chuu,  Ite  A.  Yu,  G.Y.  Yin,  J.  M.  Kahn,  and  S.  E. 
Harris,  “Hiding  Single  Photons  with  Spread  Spectrum  Technology,”  Phys.  Rev.  Lett., 
104,  223601,  (June  4,  2010). 

(2)  S.  Sensarn,  G.  Y.  Yin,  and  S.  E.  Harris,  “Generation  and  Compression  of  Chirped 
Biphotons,”  Phys.  Rev.  Lett.,  104,  253602,  (June  25,  2010). 

(3)  S.  Shwartz  and  S.E.  Harris,  “Polarization  Entangled  Photons  at  X-Ray  Energies,”  Phys. 
Rev.  Lett.  106,  080501,  2011. 

(4)  Chih-Sung  Chuu  and  S.E.  Harris,  “  Ultrabright  Backward  Wave  Biphoton  Source”,  Phys. 
Rev.  A,  83,  061803R,  2011. 

(5)  S.  Shwartz,  R.N.  Coffee,  J.M.  Feldkamp,  Y.  Feng,  J.B.  Hastings,  G.Y.  Yin,  and  S.E. 
Harris,  “X-Ray  Parametric  Down- Conversion  in  the  Langevin  Regime” ,  Phys.  Rev.  Lett., 
109,  013602,  (July  6,  2012). 

(6)  Chih-Sung  Chuu,  G.  Y.  Yin,  and  S.  E.  Harris,  “A  miniature  ultrabright  source  of  tem¬ 
porally  long,  narrowband  biphotons,”  Appl.  Phys.  Lett.,  101,  Q51108  (2012). 

(7)  T.  E.  Glover,  D.  M.  Fritz,  M.  Cammarata,  T.  K.  Allison,  Sinisa  Coh,  J.  M.  Feldkamp, 
H.  Lernke,  D.  Zhu,  Y.  Feng,  R.  N.  Coffee,  M.  Fuchs,  S.  Ghimire,  J.  Chen,  S.  Shwartz, 
D.  A.  Reis,  S.  E.  Harris  &  J.  B.  Hastings,  “X-ray  and  optical  wave  mixing,  Nature,  488, 
603,  (29  August  2012). 

(8)  S.  Shwartz,  M.  Fuchs,  J.  B.  Hastings,  Y.  Inubushi,  T.  Ishikawa,  T.  Katayama,  D.  A.  Reis, 
T.  Sato,  K.  Tono,  M.  Yabashi,  and  S.  E.  Harris, “X-Ray  Second  Harmonic  Generation  ”, 
(submitted  for  publication). 

(9)  Steven  Sensarn  “Nonlocal  Modulation  and  Dispersion  ”,  PhD  dissertation,  (May  2010). 

(10)  Chinmay  Belthangady  “Amplitude  and  Phase  Modulation  of  Single  Photons  and  Bipho¬ 
tons”,  PhD  dissertation,  (June  2010). 

The  paper  “Polarization  Entangled  Photons  at  X-Ray  Energies”  was  upgraded  by  the 
reviewers  from  a  PRA  to  Physical  Review  Letters,  and  the  paper  “Ultrabright  Backward 
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Wave  Biphoton  Source”  was  selected  by  the  editors  for  a  Physics  Viewpoint  (spotlighting 
exceptional  research). 


2.2  Presentations 

(1)  S.  Shwartz  and  S.  E.  Harris,  “Generation  of  Polarization  Entangled  Photons  at  X-Ray 
Energies  ”,  SPRC  2010  Annual  Symposium,  Stanford,  CA  (September  13,  2010)  (invited). 

(2)  S.  E.  Harris,  “Modulation  of  Photons  and  Biphotons,  Texas  A&M,  November  2010. 

(3)  Chili-Sung  Chuu,  Chinmay  Belthangady,  Ite  A.  Y,  G.Y.  Yin,  J.  M.  Kahn,  and  S.E.  Har¬ 
ris,  ’’Spread  Spectrum  Technology  with  Single  Photons  and  Slow  Light,  SPIE  Photonics 
West,  San  Francisco,  CA  (January  22-27,  2011)  [invited]. 

(4)  Chih-Sung  Chuu,  Chinmay  Belthangady,  Ite  A.  Y,  G.Y.  Yin,  J.  M.  Kahn,  and  S.E. 
Harris,  Hiding  Single  Photons  with  Spread  Spectrum  Technology,  42nd  Annual  Meeting 
of  the  Division  of  Atomic,  Molecular, and  Optical  Physics  (APS-DAMOP-2011),  Atlanta, 
GA  (June  13,  2011). 

(5)  S.  Shwartz  and  S.  E.  Harris,  ’’Polarization  Entangled  Photons  at  X-Ray  Energies”,  OSA 
Topical  Meeting  on  Nonlinear  Optics  (NLO),  Kauai,  HI  (July  17,  2011). 

(6)  Chih-Sung  Chuu  and  Steve  Harris,  ’’Ultra-bright  Backward  Wave  Biphoton  Source”, 
OSA  Topical  Meeting  on  Nonlinear  Optics  (NLO),  Kauai,  HI  (July  17,  2011). 

(7)  S.E.  Harris,  Modulation  of  Photons  and  Biphotons,  OSA  Topical  Meeting  on  Nonlinear 
Optics  (NLO),  Kauai,  HI  (July  17,  2011)  [plenary], 

(8)  S.  Shwartz,  ”X-ray  Nonlinear  Optics  in  Crystals  ”,  LCLS/SSRL  users’  meeting  SLAC, 
CA  (2011)  (Invited). 

(9)  S.  Shwartz,  R.  N.  Coffee,  J.  M.  Feldkamp,  Y.  Feng,  J.  B.  Hastings,  G.  Y.  Yin,  and 
S.  E.  Harris,  ” X-Ray  Parametric  Down-Conversion  in  The  Langevin  Regime”,  Winter 
Colloquium  on  the  Physics  of  Quantum  Electronics  (PQE),  Snowbird,  Utah  (January 
2012)  (Invited). 
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(10)  Chih-Sung  Chuu,, Generation  and  Modulation  of  Temporally  Long  Single  Photons  and 
Biphotons,  Physics  colloquium.  Physics  Department,  National  Tsing  Hua  University, 
Hsinchu,  Taiwan,  February  22,  2012. 

(11)  Chih-Sung  Chuu,,  Manipulating  the  wavepackets  of  single  photons  and  biphotons, 
Physics  colloquium,  Physics  Department,  National  Chung  Cheng  University,  Chiayi,  Tai¬ 
wan,  March  21,  2012. 

(12)  S.  Shwartz  and  S.  E.  Harris,  ”  Parametric  Down-Conversion  at  X-Ray  Wavelengths”, 
International  Workshop  on  Atomic  Physics  2012,  Dresden,  Germany  (November  2012) 
(Invited). 

(13)  Chih-Sung  Chuu,,  Shaping  single  photons  and  biphotons,  Colloquium,  Electrophysics 
Department,  National  Chia  Tung  University,  Hsinchu,  Taiwan,  November  15,  2012. 

(14)  Chih-Sung  Chuu,, Shaping  single  photons  and  biphotons,  Colloquium,  Physics  Depart¬ 
ment,  National  Taiwan  Normal  University,  Taipei,  Taiwan,  December  19,  2012. 

(15)  Chih-Sung  Chuu  Shaping  single  photons  and  biphotons,  Colloquium,  Physics  Depart¬ 
ment,  National  Cheng  Kung  University,  Tainan,  Taiwan,  January  4,  2013. 

(16)  Chih-Sung  Chuu,  Shaping  single  photons  and  biphotons,  Colloquium,  Institute  of 
Physics,  National  Chia  Tung  University,  Hsinchu,  Taiwan,  January  3,  2013. 

(17)  S.  Shwartz,  M.  Fuchs,  J.  B.  Hastings,  Yuichi  Inubushi,  T.  Ishikawa,  T.  Katayama,  D. 

A.  Reis,  T.  Sato,  K.  Tono,  M.  Yabashi,  and  S.  E.  Harris,  Second  Harmonic  Generation  at 
X-Ray  Wavelengths,  OSA  Topical  Meeting  on  Nonlinear  Optics,  Kona,  HI  (July  2013). 
(Invited) 

(18)  S.E.  Harris  “Fifty  Years  of  Parametric  Down  Conversion:  From  Microwaves  to  x-rays, 
CLEO-QELS,  San  Jose  ,  CA,  (June  2013  [plenary]). 

2.3  Patent  Disclosure 

“Ultrabright  long  biphoton  generation  with  non-linear  optical  material” ,  Sll-432 /US-PID:  18566 
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3  Abstracts  of  Published  Papers 

3.1  Hiding  Single  Photons  with  Spread  Spectrum  Technology 

Authors:  Chinmay  Belthangady,  Chih-Sung  Chuu,  Ite  A.  Yu,  G.Y.  Yin,  J.M.  Kalin,  and 
S.E.  Harris 

We  describe  a  proof-of-principal  experiment  demonstrating  the  use  of  spread  spectrum 
technology  at  the  single  photon  level.  We  show  how  single  photons  with  a  prescribed  tem¬ 
poral  shape,  in  the  presence  of  interfering  noise,  may  be  hidden  and  recovered. 

3.2  Generation  and  Compression  of  Chirped  Biphotons 

Authors:  S.  Sensarn,  G.  Y.  Yin,  and  S.  E.  Harris 

We  describe  an  experiment  demonstrating  the  radarlike  technique  of  chirp  and  compress. 
Chirped  biphotons  are  generated  using  a  quasi-phase-matched  nonlinear  crystal  where  the 
phase-matched  frequency  varies  linearly  with  position.  Sum  frequency  generation  is  used  to 
measure  the  amplitude  of  the  biphoton  wave  function.  As  compared  to  a  nonchirped  crystal, 
compression  and  an  increase  in  summing  efficiency  of  a  factor  of  5  is  observed. 

3.3  Polarization  Entangled  Photons  at  X-Ray  Energies 

Authors:  S.  Shwartz  and  S.E.  Harris 

We  show  that  polarization  entangled  photons  at  x-ray  energies  can  be  generated  via 
spontaneous  parametric  down-conversion.  Each  of  the  four  Bell  states  can  be  generated  by 
choosing  the  angle  of  incidence  and  polarization  of  the  pumping  beam. 

3.4  Ultrabright  Backward  Wave  Biphoton  Source 

Authors:  Chih-Sung  Chuu  and  S.E.  Harris 

We  calculate  the  properties  of  a  biphoton  source  based  on  resonant  backward-wave  spon¬ 
taneous  parametric  down  conversion.  We  show  that  the  biphotons  are  generated  in  a  single 
longitudinal  mode  having  a  subnatural  linewidth  and  a  Glauber  correlation  time  exceeding 
65  ns. 
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3.5  A  miniature  ultrabright  source  of  temporally  long,  narrow- 
band  biphotons 

Authors:  Chili-Sung  Chuu,  G.  Y.  Yin,  and  S.  E.  Harris, 

We  demonstrate  a  miniature  source  of  long  biphotons  utilizing  the  cluster  effect  and 
double-pass  pumping  in  a  monolithic  doubly  resonant  parametric  down-converter.  We  obtain 
a  biphoton  correlation  time  of  17.1ns  with  a  generation  rate  of  1.1  x  105  biphotons/(s  mW) 
and  an  estimated  linewidth  of  8.3  MHz. 

3.6  X-Ray  Parametric  Down- Conversion  in  the  Langevin  Regime 

Authors:  S.  Shwartz,  R.N.  Coffee,  J.M.  Feldkamp,  Y.  Feng,  J.B.  Hastings,  G.Y.  Yin,  and 
S.E.  Harris 

We  experimentally  and  theoretically  study  the  coincidence  count  rate  for  down-converted 
x-ray  photons.  Because  of  photoionization,  parametric  down-conversion  at  x-ray  wavelengths 
generally  involves  loss  and  the  theoretical  description  requires  a  Langevin  approach.  By 
working  in  a  transmission  geometry  (Laue)  rather  than  in  the  Bragg  geometry  of  previ¬ 
ous  experiments,  we  obtain  an  improvement  in  the  signal-to-noise  ratio  of  12.5,  and  find 
agreement  between  experiment  and  theory. 

3.7  X-ray  and  optical  wave  mixing 

Authors:  T.  E.  Glover,  D.  M.  Fritz,  M.  Cammarata,  T.  K.  Allis,  Sinisa  Coh,  J.  M.  Feldkamp, 
H.  Lernke,  D.Zhu,  Y.  Feng,  R.  N.  Coffee,  M.  Fuchs,  S.  Ghimire,  J.  Chen,  S.  Shwartz,  D.  A. 
Reis  ,  S.E. Harris,  and  J.  B.  Hastings 

Lightmatter  interactions  are  ubiquitous,  and  underpin  a  wide  range  of  basic  research  fields 
and  applied  technologies.  Although  optical  interactions  have  been  intensively  studied,  their 
microscopic  details  are  often  poorly  understood  and  have  so  far  not  been  directly  measurable. 
X-ray  and  optical  wave  mixing  was  proposed  nearly  half  a  century  ago  as  an  atomic-scale 
probe  of  optical  interactions  but  has  not  yet  been  observed  owing  to  a  lack  of  sufficiently 
intense  X-ray  sources.  Here  we  use  an  X-ray  laser  to  demonstrate  X-ray  and  optical  sum- 
frequency  generation.  The  underlying  nonlinearity  is  a  reciprocal-space  probe  of  the  optically 


induced  charges  and  associated  microscopic  fields  that  arise  in  an  illuminated  material.  To 
within  the  experimental  errors,  the  measured  efficiency  is  consistent  with  first-principles 
calculations  of  microscopic  optical  polarization  in  diamond.  The  ability  to  probe  optical 
interactions  on  the  atomic  scale  offers  new  opportunities  in  both  basic  and  applied  areas  of 
science 

4  X-Ray  Second  Harmonic  Generation 

Our  work  on  second  harmonic  generation  has  now  been  submitted  for  publication  .  The 
authors  of  the  submitted  publication  are:  S.  Shwartz,  M.  Fuchs,  J.  B.  Hastings,  Y.  Inubushi, 
T.  Ishikawa,  T.  Katayama,  D.  A.  Reis,  T.  Sato,  K.  Tono,  M.  Yabashi,  S.  Yudovich,  and 
S.  E.  Harris.  An  excerpt  and  summary  of  this  paper  follows: 

4.1  Introduction 

This  work  reports  experimental  evidence  for  coherent  nonlinear  x-ray  interaction  at  intensi¬ 
ties  on  order  of  10  16  W/cm2,  where  the  corresponding  electric  field  approaches  the  atomic- 
field  strength.  Using  a  1.7  Angstrom  pump  beam  we  observed  second  harmonic  generation  in 
diamond  without  permanent  damage.  The  generated  second  harmonic  is  of  order  ten  times 
the  background  radiation,  scales  quadratically  with  pump  pulse-energy,  and  is  generated  over 
a  narrow  phase-matching  condition.  Our  observation  of  a  nonlinear  effect  at  intensities  that 
are  several  orders  of  magnitude  higher  than  the  radiation  damage  threshold  in  the  visible 
regime  forms  a  basis  for  more  general  nonlinear  experiments  in  the  area  of  x-ray  nonlinear 
optics. 

It  is  now  40  years  since  Freund  and  Eisenberger  and  colleagues  described  a  theory  for  the 
nature  of  a  solid-state  dense  plasma  nonlinearity  that  is  operative  at  x-ray  wavelengths.  This 
nonlinearity  is  nonlocal  and  of  the  second  order  and  may  be  observed  in  centro-symmetric 
materials  but  requires  a  nonuniform  electron  density.  Because  of  the  small  magnitude  of  this 
nonlinearity,  to  date,  all  experiments  have  been  based  on  parametric  down  conversion  from 
the  vacuum.  The  key  factor  for  those  experiments  is  the  large  number  of  vacuum  fluctuation 
modes  at  x-ray  wavelengths.  Thus,  although  the  parametric  gain  is  very  small,  coincidences 
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of  the  generated  signal  and  idler  photons  can  be  measured  even  when  the  pump  source  is  an 
x-ray  tube.  Recently,  Glover  and  colleagues  have  demonstrated  the  nonlinear  wave  mixing  of 
x-ray  and  near-infrared  beams  using  an  x-ray  FEL.  However,  the  efficiency  in  this  depends 
on  the  intensity  of  the  infrared  laser  and  not  on  the  intensity  of  the  x-ray  laser.  On  the 
other  hand,  the  high  peak  power  an  x-ray  FEL  is  critical  for  the  observation  of  most  x-ray 
coherent  nonlinear  effects. 

Here  we  report  the  phase-matched  second  harmonic  generation  at  hard  x-ray  wavelengths. 
Using  a  1.7  A(7.3  keV)  pump  beam  from  the  SACLA  FEL  at  an  intensity  of  1016  W/cm2 
the  observed  second  harmonic  is  generated  in  a  crystal  of  diamond.  The  observed  second 
harmonic  (SH)  signal  is  about  1  photon  every  150  x-ray  laser  shots  each  containing  about  a 
billion  7.3  keV  photons  focused  to  a  micron-sized  spot  on  the  crystal.  The  second  harmonic 
beam  is  generated  in  a  narrow  angular  range  of  0.2  rnrad  full-width  at  half  maximum 
(FWHM)  and  scales  quadratically  with  pump-pulse  energy. 

The  weakness  of  the  process,  even  with  focused  x-ray  fields  on  the  order  of  the  atomic- 
field  strength,  requires  careful  attention  to  reduce  or  eliminate  potential  backgrounds.  In 
particular,  we  designed  the  experiment  to  suppress  the  off-axis  second  harmonic  in  the  pump 
beam  generated  in  the  free-electron  laser  process  itself  as  well  as  well  parasitic  scatter  of 
the  fundamental  pump.  At  hard  x-ray  wavelengths,  phase  matching  can  be  achieved  by 
using  the  periodic  nature  of  the  electron  density  in  crystals.  Phase  matching  occurs  when 
2 ku  +  G  =  k2uj  where  are  the  internal  wavevectors  of  the  pump  and  second  harmonic  beams 
respectively  and  G  is  a  reciprocal  lattice  vector.  Because  of  the  finite  dispersion,  the  angles 
between  the  beams  and  the  atomic  planes  differ  slightly  from  the  linear  Bragg  diffraction 
condition  of  the  second  harmonic.  However,  the  expected  offsets  are  small  compared  to  both 
the  divergence  of  the  focused  pump  beam  and  the  measured  width  of  the  SH  signal.  This 
presents  a  challenge  for  distinguishing  SHG  from  the  linear  Bragg  diffraction  of  the  residual 
SH  in  the  pumping  beam. 
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Figure  1:  Schematic  of  the  SHG  experimental  setup.  X-ray  pulses  from  the  SACLA  free  elec¬ 
tron  laser  are  incident  a  silicon  (111)  channel  cut  monochromator  at  10Hz.  The  monochroma¬ 
tor  selects  a  7.3  keV  pump-beam  with  leV  bandwidth  while  rejecting  a  substantial  fraction 
of  the  second  harmonics.  A  pair  of  Kirkpatrick  Baez  mirrors  focuses  the  beam  to  about  1.5 
microns  while  reducing  the  third  harmonic  contamination.  The  nonlinear  medium  is  a  (111) 
cut  diamond  crystal  placed  in  the  focus  and  set  for  phase  matching  as  indicated  in  the  phase 
matching  diagram  (bottom).  The  phase  matching  is  achieved  by  using  the  (220)  atomic 
planes;  k ^  and  k2U  are  the  wave  vectors  of  the  pump  and  the  generated  second  harmonic 
respectively.  G  is  the  reciprocal  lattice  vector  orthogonal  to  the  (220)  atomic  planes.  We 
control  the  intensity  of  the  pump  by  using  a  set  of  A1  filters  before  the  sample. 
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4.2  Experiment 


The  experiment  was  performed  at  BL3  SPring-8  Angstrom  Compact  free  electron  laser 
(SACLA)  in  Hyogo,  Japan.  The  experimental  setup  is  shown  schematically  in  Fig  1. 

In  order  to  observe  the  SHG  generated  in  the  crystal  extreme  care  was  required  to  sup¬ 
press  backgrounds  from  residual  second  harmonic  contribution  in  the  beam  and  parasitic 
scattering  of  the  fundamental  into  the  detector  as  well  as  orient  the  sample  and  detector.  A 
Si  (111)  double-crystal  monochromator  is  used  to  select  a  narrow  (leV)  bandwidth  pump- 
beam  centered  at  7.3  keV  photon  energy  from  the  FEL.  While  even  harmonics  are  forbidden 
by  symmetry  from  being  radiated  on  axis,  there  is  a  finite  SH  generated  slightly  off-axis 
due  to  the  highly  relativistic  motion  of  the  electrons  in  the  undulator.  The  monochromator 
suppresses  the  second  harmonic  by  more  than  four  orders  of  magnitude  because  the  second 
order  Si  (222)  Bragg  peak  is  very  nearly  forbidden.  After  monochromatization  the  beam 
is  focused  by  a  pair  of  grazing  incidence  mirrors  in  the  Kirkpatrick-Baez  (KB)  geometry  to 
about  1.5  m  FWHM  spot  at  the  sample  (assuming  a  Gaussian  beam,  the  Rayleigh  range  is 
30  mm  and  the  divergence  is  41  rad).  A  set  of  thin  A1  filters  before  the  KB  mirrors  are  used 
to  attenuate  the  fundamental  while  having  negligible  effect  on  the  residual  harmonic.  The 
pulse  energy  after  the  filters  is  measured  on  a  shot-by-shot  basis.  The  sample  is  a  0.48  mm 
thick  (111)  cut  diamond  single  crystal  and  is  placed  in  the  focus  of  the  x-ray  beam.  The 
sample  is  oriented  near  the  diffraction  condition  for  the  (220)  planes  in  the  symmetric  Laue 
geometry  with  the  polarization  of  the  pump  in  the  scattering  plane.  In  this  geometry  both 
diffracted  residual  second  harmonic  and  the  generated  second  harmonic  transmit  through 
the  crystal.  Scattered  photons  are  detected  by  a  YAP:Ce  scintillation  detector  with  energy 
resolution  of  30.  We  choose  to  phase-match  with  the  (220)  reciprocal  lattice  vector.  Impor¬ 
tantly,  this  corresponds  to  the  (110)  planes  for  the  fundamental  which  is  strictly  forbidden 
in  the  face-center  cubic  structure,  and  thus  elastic  scattering  of  the  pump  in  the  direction  of 
the  second  harmonic  generation  is  very  highly  suppressed.  Nonetheless,  as  we  expect  average 
count  rates  of  much  less  than  a  single  second  harmonic  photon  per  FEL  pulse,  additional 
suppression  is  necessary.  Note  that  the  detector  cannot  distinguish  a  single  SH  photon  from 
two  fundamental  photons  on  any  given  pulse.  Moreover,  the  probability  of  measuring  two 
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photons  per  pulse  also  scales  quadratically  with  incident  pulse  energy  (as  is  the  expected 
SHG  signal),  so  the  suppression  of  stray  scattered  fundamental  photons  is  critical.  We  use 
slits  and  a  set  of  A1  foils  between  the  sample  and  the  detector  to  restrict  the  solid  angle 
of  the  detector  and  filter  out  background  photons  at  the  fundamental  wavelength.  The  A1 
filters  attenuate  the  background  at  the  fundamental  wavelength  by  a  factor  of  106,  while  the 
SHG  beam  is  attenuated  by  a  factor  of  6.75. 

4.3  Results 

Since  the  maximum  intensity  in  our  experiment  was  about  1016  W/cm2,  we  first  checked  for 
damage  by  measuring  the  elastic  scattering  of  the  (220)  peak  at  7.3  KeV  at  the  maximum 
pulse  energy  and  observed  no  changes  in  the  peak  reflectivity  or  in  the  rocking  curve.  We 
conclude  that  the  crystal  structure  did  not  change  during  the  exposure  to  the  XFEL  beam. 

A  histogram  of  detector  counts  as  a  function  of  photon  energy  is  shown  in  Fig.  2a  at 
full  intensity  (no  filter  before  the  diamond  crystal)  and  at  the  peak  of  the  phase  matching 
condition  for  24000  shots.  The  energy  bin  width  is  1.39  keV.  The  data  of  Fig  2a  consists 
primarily  of  two  peaks  corresponding  to  photon  energies  near  the  fundamental  and  SH.  There 
is  a  small  possibility  of  detected  third  harmonic  photons  as  well.  If  we  apply  upper  and  lower- 
level  thresholds  for  defining  a  SH  photon  as  10. 2-19  KeV  (+/-  30%  of  14.6  keV),  we  measure 
153  SH  photons  in  24000  pulses  (  0.00640.0005/pulse),  compared  to  35  fundamental  photons 
in  5.19.5  keV  (0.00150.0002/pulse).  This  level  is  low  enough  that  contamination  from  the 
wings  of  the  fundamental  (or  third  harmonic)  is  small,  and  the  probability  of  two  fundamental 
photons  in  a  single  pulse  is  negligible. 

In  order  to  determine  whether  the  second  harmonic  photons  are  generated  in  the  crystal  or 
due  to  elastic  scattering,  We  insert  a  0.0025  mm  A1  filter  before  the  diamond  crystal  in  order 
to  determine  whether  the  SH  photons  are  generated  in  the  crystal  or  due  to  elastic  scattering 
from  the  residual  harmonic  content  in  the  incident  beam.  The  absorption  coefficients  for  A1 
at  the  fundamental  and  SH  are  168.0  cm-1  and  21.8  crn-1  corresponding  to  transmission 
through  the  pre-filter  of  0.43  and  0.9  respectively.  The  resulting  histogram  is  plotted  in  Fig. 
2b.  In  this  case  we  count  only  7  second  harmonic  photons  in  6000  pulses,  207  %  of  what 
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one  would  expect  if  the  signal  was  solely  from  undulator  SH.  This  indicates  that  most  of 
the  signal  at  the  highest  pump  pulse-energy  is  SH  generated  in  the  crystal.  In  fact,  if  the 
SH  signal  scales  quadratically  with  pulse  energy,  than  the  signal  upon  inserting  the  A1  filter 
should  be  reduced  by  a  factor  of  slightly  more  than  five  to  18.5%. 

To  better  quantify  the  background  we  repeat  the  measurement  with  0.025  mm,  0.1  mm, 
and  0.2  mm  thick  A1  filters  and  with  no  filter.  Fig.  3a  shows  the  SH  count  rate  as  a  function 
of  the  average  of  the  pump  pulse-energy  when  the  crystal  is  detuned  by  -57  rad  from  the 
peak  of  the  phase  matching  condition.  Fig  3b  shows  the  rocking  curve  of  the  SHG  process, 
namely,  the  SH  count  rate  as  a  function  of  the  angular  deviation  of  the  crystal  from  the 
phase  matching  angle  taken  without  filters.  Consistent  with  the  data  on  the  peak  of  the 
phase  matching  condition,  the  SH  signal  at  -57  rad  (Fig  3a)  upon  inserting  the  0.025  mm 
pre-filter  is  reduced  to  188  %  of  the  value  without  the  filter.  It  is  clear  from  Fig.  3a  that  the 
count  rates  at  the  two  lowest  pulse  energies  are  equal  within  the  statistical  error.  We  use 
this  to  estimate  the  residual  SH  assuming  that  it  dominates  the  measurement  at  the  lowest 
pulse-energy.  After  correcting  for  the  filter  transmission,  we  find  that  the  residual  harmonic 
corresponds  to  0.00540.0030  counts/pulse.  Thus  the  SHG  from  the  crystal  is  about  a  factor 
of  ten  higher  than  the  background  at  the  highest  intensity.  Assuming  that  the  background 
is  negligible  at  all  angles,  the  width  of  the  curve  represents  the  width  over  which  the  second 
harmonic  is  phase  matched.  We  plot  a  Gaussian  fit  and  find  the  FWHM  of  the  rocking  curve 
is  about  180  p  rad. 

In  summary,  we  report  experimental  evidence  for  phase  matched  SHG  at  hard  x-ray 
wavelengths.  The  magnitude  of  the  generated  SH  energy  scales  as  the  square  of  the  pump 
pulse  energy,  and  its  magnitude  peaks  at  the  calculated  phase  matching  angle.  The  largest 
observed  SHG  efficiency  is  5.8  x  ICG11,  with  an  observed  count  rate  that  is  more  than  ten 
times  above  the  background  noise.  These  results  advance  our  understanding  of  coherent 
x-ray  matter  interactions  at  high  intensity  and  form  the  basis  for  more  general  experiments 
on  x-ray  nonlinear  processes.  For  example,  future  experiments  might  explore  x-ray  phase 
conjugation,  squeezing,  and  the  generation  of  polarization  entangled  x-ray  photons. 
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Figure  2:  Photon-energy  histograms  of  the  measured  signal  at  the  peak  of  the  phase-matching 
condition.  The  pump-photon-energy  is  7.3  keV  and  the  SHG  photon-energy  is  14.6  keV.  The 
energy  resolution  of  the  detector  is  30  %.  The  data  in  part  (a)  are  taken  with  no  attenuation. 
The  data  in  part  (b)  are  taken  with  a  0.025  mm  A1  filter  before  the  diamond  crystal  (the 
transmission  through  the  filter  is  0.43  and  0.9  at  7.3  keV  and  at  14.6  keV  respectively). 
In  Fig.  2b  the  second  harmonic  count  rate  is  20  percent  of  what  would  be  expected  if  it 
resulted  from  the  undulator  alone.  Note  that  the  number  counts  per  pulse  is  small,  and  the 
second  harmonic  counts  are  larger  than  the  counts  at  the  fundamental  frequency;  thus,  the 
contribution  from  two  simultaneous  fundamental  photons  is  negligible. 
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Figure  3:  Scaling  of  second  harmonic  rate  with  average  pump  pulse-energy  and  deviation 
angle,  (a)  Detected  second  harmonic  signal  as  a  function  of  input  pulse-energy  and  measured 
at  angular  deviation  of  -57  /i  rad  from  the  SHG  phase-matching  angle.  The  pulse-energy 
is  varied  by  inserting  thin  A1  filters  in  the  beam  before  the  diamond  crystal.  Each  of  the 
measured  data  points  represents  an  average  over  6000  pulses.  The  dashed  blue  curve  is 
calculated  from  the  theory  and  scaled  vertically  by  a  factor  of  1.48.  The  second  harmonic 
contribution  at  the  higher  pulse  energies  (two  points)  is  dominated  by  SHG,  while  at  the 
lower  pulse  energies  the  residual  second  harmonic  contamination  dominates,  (b)  Normalized 
second  harmonic  signal  as  a  function  of  the  angular  deviation  of  the  diamond  crystal  from 
the  phase  matching  angle.  The  data  in  part  (b)  are  obtained  with  no  attenuation  before  the 
diamond  crystal.  The  dashed  curve  is  a  Gaussian  fit  to  the  experimental  data.  The  FWHM 
of  the  rocking  curve  is  180  /i  rad.  The  peak  point  represents  the  average  signal  measured 
over  2400  pulses.  The  other  points  represent  the  average  signal  measured  over  6000  pulses. 
In  both  plots,  the  vertical  error  bars  indicate  the  counting  statistics. 
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5  Temporal  and  Spatial  Correlation  of  x-ray  biphotons 


It  is  the  intent  of  this  Section  to  derive  and  summarize  the  formulli  for  the  temporal  and 
spatial  correlations  that  are  inherent  to  the  spontaneous  downconversion  of  an  incident  near- 
monochromatic  x-ray  pumping  beam.  When  combined  with  further  numerical  work,  it  is 
likely,  but  not  certain  that  we  will  submit  the  work  of  this  section  for  publication.  This 
section  was  written  by  S.E.  Harris. 

We  begin  by  developing  the  coupled  equations  for  this  inherently  non-collinear  problem. 
To  start,  we  separate  the  driving  term  of  the  wave  equation  into  a  polarization  V  that  is 
linear  in  the  driving  e-field  and  a  driving  current  density  J  that  is  quadratic  in  the  driving 
e-fields.  With  Vj _  =  ax-£:  +  ay-^,  the  wave  equation  in  free  space  is 


d2£  2  1  d2£  _  d2V  dj 

+  V±£  ~  W  +  /i0DP 


(5.1) 


With  r  =  {x,y},  and  q  =  {kx,ky},  we  Fourier  Transform  Eq.  (1)  so  that  £(z,x,y,t)  —> 
E(z,  kx,  ky,uj),  V(z,x,y,t)  — >  P(z,  kx,  ky,u> )  and  £f(z,x,y,t)  —$■  J{z ,  kx,  ky,  u).  For  example, 


£(z,r,t)  = 


E(z,  q,  lo)  exp[— i(ujt  —  q  •  r)]dqdu; 


(5.2) 


'  —oo  J  —oo 


We  substitute  Eq.(2)  into  Eq.(l)  and  write  the  linear  portion  of  the  frequency  domain  dipole 
moment  in  terms  of  the  refractive  index  P  (u)  =  £o[n(u})2  —  l]E(u}).  With/c2(a;)  =  ui2n2(ui) /c2, 
the  Fourier  transform  of  Eq.(l)  is 

d2  E 

q  2 —  T  ky)E  +  k2(uj)E  =  — .  (5-3) 

Eq.(3)  is  exact  and  contains  dispersion  to  all  orders. 

It  is  often  the  case  that  the  electromagnetic  fields  vary  slowly  in  the  z-direction  as  com¬ 
pared  to  the  x-ray  wavelength.  It  is  then  useful  to  the  change  variable  from 


E  =  E(z)  exp (ikzz)  (5.4) 

J  =  J(z)  exp (ikzz) 

With  k2  +  ky2  +  k2  =  k2(uo)  =  uj2n(u)2  /  c2 ,  and  assuming  that  d2E/dz 2  «  kzdE/dz  we 
obtain  the  equation  for  the  slowly  varying  Fourier  components  of  the  electromagnetic  fields. 


17 


With  the  impedance  r](u)  =  \J (/ioAo)^(^),  and  an  angle  of  propagation  9Z  with  regard  to 
the  z  axis,  we  have 

cos  (9z)dE/dz  =  — — (5.5) 

In  the  following  we  will  assume  that  the  signal  and  idler  fields  are  described  by  the  proto- 
type  plasma  nonlinearity  as  first  described  by  Levine  and  Eisenberger  [] .  This  nonlinearity  is 
dependent  on  whether  the  various  fields  are  polarized  in,  or  out,  of  the  plane  of  incidence  [], 

and  generally  leads  to  the  generation  of  polarization  entangled  photons  [] .  This  nonlinearity 

,  —  ~  ~  ^  ~  ~  ~  ^ 

satisfies  detailed  balance,  so  that  with  ((us,Ui,Up)  =  (,  Js  =  usQEpEi  ,  A  =  c Oi(EpEs  ,  and 
Jp  =  —ujp(EsEi.  Following  the  transformation  of  Eq.(4),  the  nonlinearity  as  it  appears  in 
the  first  of  the  coupled  S.V.E.A.  equations  is  Js  =  u js(EpEi*  exp(iAk~z)  where 

Akz  kpZ  ( ksz  T  kiz)  (5.6) 


The  coupled  S.V.E.A.  equations  for  the  signal,  idler,  and  pump  fields  are 


Zp  (]  71 

cos  (Os)-rp-  +  asEs  = - ~^—EpE*  exp(iAkzz) 

dE'  (jj'Ti'C 

cos (9i)-rpr  +  aiEi  = - Tp—EpEg*  exp(iAkzz) 

ta  \  ®EP  i  cl  UpllpC, 

cos(vp)^T  +  aPEP  = 


dz 


EsEi  exp(—iAkzz) 


(5.7) 


With  the  loss  at  all  frequencies  taken  as  zero,  and  with  9  as  the  angle  with  the  z-axis  these 
equations  satisfy  the  photon  conservation  condition 


cos(d 
—  cos[9p 


1  d  (\E 


los  dz  \  2rj, 


=  cos  (9i) 


1  d  ( \Ei 


u>i  dz  \  2 rj, 


1 2  \  1 


1  d 

fWl 

Up  dz 

K2Vp  J . 

(5.8) 


5.1  Coupled  Operator  Equations 


We  work  in  the  Heisenberg-Langevin  picture  and  assume  that  the  pumping  field  is  indepen¬ 
dent  of  z,  and  convert  Eq.(7)  from  E-field  to  operator  notation  So  that  the  right  hand 
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side  of  each  of  Eqs.(7)  is  the  same,  we  take  Et  ~  y/oJiVi/ cos (6i)  The  coupled  operator 
equations  for  the  Fourier  components  are  then 


+  —£z~bs  =  Kb\  esxp[iAkzz]  +  J  —^rfs, 

OZ  COS  us  V  COS  o s 

db]  on  ,+  r  ■  a  t  t  .  /  2 ai  t 


~a 1  H - =  K*bsexp[-iAkzz\  + 


dz  cos  6 


cos  9, 


(5.9) 


Here,  with  r  =  1/y/ cos (9S)  cos (6**)  cos (9P),  the  coupling  coefficient  is  k  —  y  js0Ji0Jpr]srjir]p. 
Before  proceeding  with  the  solution,  we  transform  Eqs  (9)  to  constant  coefficients  by  bs  = 
as  exp(i Akzz/2)  and  b]  =  a\  exp(— iAkzz/2).  We  then  have 


das 

dz 

da\ 

dz 


+ 


+ 


Ole 


cos  9f 

Ctj 

cos  9, 


+  i 


.AK 
2 
A  k. 


—  i- 


as  =  K,a\  + 


a]  =  n*as  + 


2  as 
cos  9^ 

2  Oii 

cos  9, 


-fs, 


■ f l 


(5.10) 


The  operators  a(z,  q,  oj)  are  the  coarse  grained  annihilation  operators  with  commutator 


[ai(zi,qi,£Ui),a[.(z2,q2,W2)]  =  p-^(zi  ~~  ^)<5(qi  -  q2)<5(^i  -  u>2).  (5.11) 

They  are  normalized  so  that  the  signal  and  idler  count  rates  per  area  are  (aj(r,  t)as( r,  t))  and 
(a\ (r,  t)a,i{ r,  t))  respectively.  The  fs(z,  q,  oj)  and  // (z,  q,  oj)  are  the  Langevin  noise  operators, 
and  satisfy 


fj(z,q,cj),fl(z',q',cj') 


?^~A3^’kS(z  ~  ^)<Kq  -  q')<*(w  -  W).  (5.12) 


We  have  shown  numerically  that  the  commutator,  Eq.(ll)  is  conserved  for  all  z.  If  the 
loss  at  both  the  signal  and  idler  frequencies  is  zero,  then  the  Langevin  terms  may  be  taken 
as  zero,  while  still  preserving  the  commutator.  But  when  the  loss  is  not  zero  these  terms 
are  essential  for  obtaining  the  correct  generation  rates  at  the  signal  and  idler,  as  well  as  the 
coincidence  count  rate. 
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5.2  Solution 


With  L  as  the  crystal  length  we  write  the  solution  of  Eq.  (10)  in  terms  of  a  matrix  M  and 
the  exponential  of  this  matrix  exp  f—ML^j 


M  = 


and 

The  solution  of  Eq.  (10)  is 
flt[n  =  exp  (-ML 

L  aAL)  J  v 

Expanding 


as  +  iAkz/2  k 

k  Oii  —  iAkz/2 


A  B 
C  D 


exp  —ML  = 


as{  0) 

4(o) 


+ 


exp  (-M(s  -  L) 


(5.13) 


fs(s) 

fils) 


ds  (5.14) 


where 


ds(L)  —  ^4ns(0)  +  Bat  (0)  +  /  E(s)fs{s)ds+ 

Jo 

rL 

ft/ 


F(s)fJ  (s)ds 


Jo 

rL 


a\(L)  —  Cas(0)  +  Da\(  0)  +  /  G(s)  f s(s)ds-\- 

J  o 

rL 

ft/ 


(s)ds 


'0 


F(s )  =  \J2a. 
G(s)  =  \[2at 
H(s)  =  \/2a. 


1  1 

! 

i _ i 

'  1  ' 
0 

M(s-L)  1  • 

'  0  ' 
1 

i  i 

'oT 

i 

'  1  " 
0 

M(s-L)  1  • 

'  0  ' 
1 

(5.15) 


(5.16) 
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5.3  Signal  and  idler  count  rates 


The  quantities  of  interest,  that  is  the  generation  rates,  coincidence  rates,  and  the  temporal 
and  spatial  correlation  functions  may  all  be  written  in  terms  of  the  frequency  domain  func¬ 
tions  that  are  defined  in  Eqs(15)  and  (16).  The  rate  Rs  [photons/(sec  area)]  of  generation 
of  signal  photons  is  given  by  Rs  =  (aj (r,  t)as(r,  t)).  We  Fourier  transform  from  t  — *  co, 
x  — >  kx,  and  y  — >  ky,  and  use  Eq.(15)  and  Eq.(16),  as  well  as  the  commutators  of  Eq.(ll) 
and  Eq.(12).  With  some  algebra  one  obtains 


(4(f,r)as(f,r)) 


(2tt) 


\B\2duidq  + 


uj  J  q 


(2  Jr) 


\F(s)\2dudqds  (5.17) 


id  J  q  J  0 


Similarly,  the  generation  rate  Rj  [photons/(  sec  area)]  at  the  idler  wavelength  is 


<«!(*.  r)a,(t,T)) 


(2*0 


\C\2dudq  + 


U)  J  q 


(27i -y 


\G(s)\2dudqds  (5.18) 


u  J  q  J  0 


Because  of  the  loss,  the  generation  rate  at  the  signal,  may  or  may  not,  be  equal  to  that  at 
the  idler  wavelength.  We  assume  a  plane  wave  pump  beam  with  a  projected  area  on  the 
x-y  plane  Ap;  so  that  the  total  count  rate  at  either  the  signal  or  the  idler  is  obtained  by 
multiplying  by  Ap.  But  note,  that  since  the  square  of  the  coupling  coefficient  in  Eq.  (10)  is 
inversely  proportional  to  Ap,  the  total  generated  count  rate  is  independent  of  this  area. 

From  their  definitions  in  Eq.(15)  we  see  that  when  the  loss  at  the  signal  frequency  is  equal 
to  zero,  the  quantity  \A\2  =  \A{uj,  kx,  ky)\2  is  the  power  gain  Gs  at  the  signal  wavelength. 
The  standard  quantum  limit  requires  that  the  count  rate  at  the  signal,  i.e.,  the  noise  per 
temporal  and  angular  bandwidth  at  the  signal  frequency,  be  equal  to  Gs(co)  —  1.  With  the 
noise  given  by  Eq.(17),  we  require  that  when  as  =  0, 

rL 


i^r 


l  =  I  By  + 


|F(s)|2ds 


(5.19) 


Similarly,  when  cq  =  0, 


\D\ 


l  =  |Cf+/  |G(s)|2ds 


(5.20) 


We  have  numerically  veriEed  that  the  quantities  A . H  satisfy  Eqs.(19)  and  (20). 
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5.4  Correlations  and  coincidence  counts 


Assuming  two  ideal  detectors  each  with  an  area  A  that  is  larger  than  that  of  the  pumping 
beam,  then  if  the  first  detector  clicks  at  time  t\  and  position  ri  =  {cci,  2/1},  the  probability 
of  the  second  detector  clicking  at  time  t2  =  t\  +  t  and  position  r2  =  ri  +  /1,  is  given  by  the 
Glauber  correlation  function  G ^ 


G,2,(r,fi)  =  {o*((2,r2)a?(i1,r,)ai(«1,r1)o,(f2,  r2)) 


(5.21) 


To  evaluate  G^(r,  /i),  we  Fourier  transform,  use  Wick’s  theorem,  Eqs.  (15)  and  (16)  and 
the  commutators  of  Eqs  (11)  and  (12).  With  q  =  {kx,kyj,  it  is  useful  to  form  the  Fourier 
domain  function 

0(<u,q)  =  B*D+  [  F*(s)H(s)ds  (5.22) 

Jo 


With  considerable  algebra  we  obtain 


G2(t)  = 


2vr 


4>(ui,  q)  exp(iu)T)duj 


dq  +  RsRi 


(5.23) 


The  background  (or  accidental  coincidence  count  rate  RSRL  is  the  result  of  signal  and  idler 
counts  that  accidentally  arrive  within  the  same  temporal  window.  Most  often  it  is  small  as 
compared  to  the  first  term.  The  coincidence  count  rate  per  area  [photons/(sec  area)]  is 

3 


I  G(t)<It  =  J  J  \<l,(u,q)\2dudq 


(5.24) 


We  may  also  evaluate  the  spatial  correlation  of  the  incident  photons.  For  example  if  x2  = 
x\  +  fi,  then  with  cj)  =  (p(uj,  kx,  ky), 


G2(/i)  = 


2tt 


UJ  J  ky 


(f)  exp(—ikxx)dkz 


dkydoj 


(5.25) 


As  a  result  of  the  much  wider  angular  aperture  in  the  out  of  plane  dimension  y,  the  correlation 
distance  in  this  direction  will  be  much  narrower  than  in  the  in-plane  dimension. 
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5.5  A kz  in  the  Laue  geometry 


In  the  Laue  geometry  as  described  here,  the  input  and  output  planes  of  the  diamond  crystal 
are  normal  to  the  z-axis  and  the  monochromatic  pump  is  in  the  x-z  plane.  We  phase  match 
using  the  (220)  atomic  planes  with  the  reciprocal  lattice  k  vector  G  in  the  direction  of  the 
x-axis.  Denoting  the  magnitude  of  the  k-vectors  of  the  signal,  idler  and  pump  as  ks,  hi, 
and  kp ,  the  magnitude  of  the  Bragg  angle  0B  is  arcsin  (G/2kp).  The  monochromatic  pump 
propagates  in  the  x-z  plane  at  angle  9P  =  —  [9B  +  A),  where  A  is  most  often  less  than  a  few 
mr.  The  phase  matching  angles  of  the  signal  6S  and  the  idler  9^  are  also  in  the  x-z  plane  and 
are  determined  by  the  solution  of  ks  +  ki  =  kp  +  G.  For  example,  for  a  pump  at  18  Kev  the 
magnitude  of  the  Bragg  angle  is  .276  radians.  For  A  =  1  mr,  9p  =  —.277  radians,  9S  =  .244 
and  6*j=.308  radians,  respectively. 

With  the  phase  matching  angles  of  the  central  rays  determined  and  denoted  by  9S  and  9i, 
we  define  kpzO  =  kpcos(9p),  ksx 0  =  kssin(9s )  and  kix 0  =  ki  sin(6>i) ;  we  take  ksy o  =  kiy o  =  0 
and  require  that  the  sum  of  the  tangential  k-vectors  in  the  x-y  plane  is  zero.  (For  example, 
kpx  +  ksx  +  kix  +  G  =  0.)  Noting  that  kpz o  =  ksz 0  +  kiZ o,  the  k- vector  mismatch  in  the 
z-direction  is  then 

A kz  =  kpz o  -  \Jk2s~  k2sx  -  k1y  -  \J kf  -  ( ksx0  +  kix0  -  ksxf  -  k2sy  (5.26) 

When  ksx  =  ksx0,  and  klx  =  kix 0,  then  A kz  =  0.  With  the  group  velocity  at  the  signal  and 
idler  defined  as  Vs  and  V),  then  to  lowest  order  in  Scos  =  ojs  —  u^o,  and  5ksx  =  ksx  —  ksxo,  the 
quantity  A k~  expands  to 

A kz  =  [sec(9i)/Vi  -  sec(6>s)/I/s]^a;s  +  [tan(0s)  -  ta n(9i)]6ksx  +  kpscc^  scc^‘^  ksy 2  (5.27) 


6  Appendix:  Cavity-Enhanced  SPDC 

This  appendix  describes  the  theory  associated  with  the  experimental  paper  “A  miniature 
ultrabright  source  of  temporally  long,  narrowband  biphotons” ,  and  was  written  by  Chili-Sung 
Chuu. 
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6.1  Nonlinear  Polarizations 


The  real  generated  polarizations  for  the  signal  and  idler  fields  are 

Ps{z,t)  =  e0x{2)[Ep(z,t)  +  Ei(z,t)}2 

Pi(z,t)  =  e0x(2)  [Ep(z,  t)  +  Es(z,  t)}2  (6.1) 

where  Ep(z,t),  Es(z,t),  and  Ei(z,t )  are  also  real  numbers. 

We  write  the  traveling-wave  pump  field  as 

Ep(z,  t )  =  ’  -  exp [i(kpz  —  c opt)]  +  c.c.  (6.2) 

where  Ep(z,t)  is  slowly-varying  envelope,  and  approximate  the  standing-wave  signal  and 
idler  fields  by 

~  ,  .  EJt)  .  qirz  .  .  . 

Es{z,t)  =  — - —  sm——  exp(-Mns0t)  +  c.c. 

Z  Lj 

iEs{t)  ,.qnz  .qnz 

= - - —  [exp(*— )  -  exp(-2— )J  exp(-^o^)  +  c.c. 

-  .  EJt)  .  rirz  .  .  . 

Ei(z,t)  =  — —  sm—  exp(— 2uy0t)  +  c.c. 

=  -lEl^  [exp(i^)  -  exp(-i^)]  exp(-ia;iot)  +  c.c.  (6.3) 

where  q  and  r  are  integers  corresponding  to  the  cold  cavity  resonant  frequencies  Qq  = 
qnc/n^L  and  h2r  =  mc/nf^  L,  respectively,  with  n ^  and  nf^  being  the  group  indices 
of  the  signal  and  idler  fields  and  L  being  the  cavity  length.  uso  =  ^(cjp  +  —  flr)  and 

—  flq  +  flr)  are  the  central  frequencies  of  the  signal  and  idler  fields,  which  are 
near  Qq  and  h2r,  respectively,  and  satisfy  the  energy  conservation,  ujp  =  uso  +  UiQ. 

Substituting  the  expressions  of  Es(z,t ),  Ei(z,t )  and  Ep(z,t )  into  Ps(z,t )  and  Pi(z,t),  we 
then  obtain 
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Ps(z,t) 

Pi(z,t) 


eoX 


(2) 


eoX 


(2) 


Ep(z,t )  ?£)(£)  r  ,.rnz.  ,  r'KZ.,  ,  .  . 

- - - exp [i(kPz  -  ujpt)\ - - —  [exp(z— )  -  exp(-*— )J  exp(-zu;io£)  +  c.c. 

exp [i(kpz  -  upt)\  -  — ^  [exp(2^)  -  exp(-i^)]  exp(-itvs0t)  +  c.(c6.^) 


We  keep  the  terms  in  Ps(z,  t )  and  Pi(z ,  t )  that  have  frequencies  of  cuso  and  (Ujo,  respectively, 
with  the  signal  and  idler  fields  propogating  in  the  same  direction  as  the  pump, 


Pa{z,t) 

Pi(z,t) 


ie  oX 


(2) 


ieoX 


(2) 


E*(t)Ep(z,t )  exp [i(kp 
E*(t)Ep(z,t )  exp [i{kp 


TIT 

—)z  -  icOsot]  +  C.c. 
Ju 

qir  v  .  n 
—  JZ  —  *CUjotJ  +  C.c. 

1j 


(6.5) 


Since  only  the  components  that  are  nonorthogonal  to  sm(qTrz/ L)  or  sin(r7rz/L)  interact 
with  the  signal  or  idler  fields,  we  project  the  polarizations  against  the  corresponding  spatial 
modes, 


Pa,q(t)  =  jQ  Ps(z,t)sUl(<P^)dz 

s,  2  *  TIT  Z 

Pi,r(t )  =  -  Pj(z,£)sin(— — )dz  (6.6) 

P  Jo  L 


For  Ep(z,t )  =  Ep(t),  the  projected  polarization  at  the  signal  frequency  is 
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Ps,q(t)  =  Y  f  — y —  E* ( t)Ep(t )  exp [i(kp  -  ^-)z  -  iujs0t\  sin(^)d 2  +  c.c. 


LJ0  4 


,(2) 


....  ^  /  "  ™v,  exp(i^)  -  exp(-i^) 


2L 


E*(t)Ep(t )  exp(-iu;soO  /  exp[2(fcp - -)z] 


2z 


dz  +  c.c. 


eoX(2)  f  exp[i(/cp  -  ^  f  )L]  -  1  exp[i(/cp  -  ^  -  ^)L\  -  1 


4  L  i(kp-rf  +  f)L  l{kp-r,_qf)L 

e0/\(2)  f  exp [i(kp  -  rf  -  f)L\  -  1 


E*(t)Ep(t )  exp(-icusot)  +  c.c. 
-E* ( t)Ep(t )  exp(— iu;so£)  +  c.c.  (6.7) 


4  l  -rf-f)L 

where,  in  the  last  equality,  we  only  keep  the  term  in  which  the  signal  propagates  in  the  same 
direction  as  the  pump.  Defining 


a  1,  7  ttt  qn 

Ak  =  k*---~ 


(6.8) 


the  projected  polarization  at  the  signal  frequency  can  be  written  as 


PS,q(t)  =  i 


.eoX 


(2) 


4 

Ps,q(t) 


exp(zAZcL)  —  1 


A  kL 
exp(—iu>s0t)  +  c.c. 


E*(t)Ep(t)  exp(—iu}s0t)  +  c.c. 


(6.9) 


where  the  slowly-varying  envelope 


Ps,q(t)  — 


eoX 


(2) 


exp(iAkL/2 )  sinc(AA;L/2)  E*(t)Ep(t) 


(6.10) 


Similarly,  the  projected  polarization  at  the  idler  frequency  is 


P%,r(t)  = 


COX 


(2) 


4 

P%,r{t) 


exp(?A/cL/2)  sinc(A/cL/2)  E*(t)Ep(t)  exp(— iuj^t)  +  c.c. 
exp(— icoiot)  +  c.c.  (6-11) 


where  the  slowly-varying  envelope 


Pi,r(t)  = 


eoX 


(2) 


exp(zA/cL/2)  sinc(A/cL/2)  E*(t)Ep(t) 


(6.12) 
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6.2  SVEA  (Lamb)  Equations 

Defining  Es(t)  and  Ei(t)  as 


Es(z,t)  =  Es(t)  sin(^) 

^  ^  VTT  7 

Ei(z,t )  =  Ei(t)  sin(— ) 

the  second-order  differential  equation  for  Zys(f)  and  Ps,q(t)  is  given  by 


d2Es(t )  nq  dEq{t) 

dt2  Qs  dt 


+  tfqEs(t) 


1  d2PS)q(t ) 

eQes  dt 2 


where  Qs  =  eoesTlq/as  is  the  cavity  quality  factor. 


Substituting, 


(6.13) 


(6.14) 


P,(t) 

E,(t) 


exp(—iujs0t)  +  c.c. 
exp(—  ius0t)  +  c.c. 


(6.15) 


into  the  second-order  differential  equation  and  picking  out  the  e  lUs0t  terms,  we  have 


d2Es(t ) 
dt 2 


Tla  .  dEJt )  ,  o  9  .  fL  .  .  . 

~Q  2zo;so)  w so)Es(t ) 


1  g2P^g(t) 
e0es  dt2 


9.  3Ps,g 

2*cu  - 


•'sCT 


(0 


dt 


(*)] 


(6.16) 


The  slowly-varying  envelope  approximations  assumes  that 


Furthermore,  since  (nso  ~  tlq, 


<92  c> 

a?  <<:  <<:  w 


2 

sO 


(6.17) 
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q  hJ so  4~  ^bo)(^2q  ^ho)  ~  2cUso(f2g  hJ sO )  (6.18) 

We  also  assume  that  the  cavity  decay  rate  is  small  compared  to  the  resonant  frequency  of 
the  cavity 

rs  =  ^  <  nq  (6.19) 


With  these  approximations,  we  obtain  the  complex  SVEA  equation 


dEs(t) 

dt 


k-bo)  + 


2Qs 


Similary,  for  Ei(t)  and  PltT,  we  have 


Ea(t) 


ic^so 


Ps,q{t) 


dEi(t) 
dt  + 


<^io)  + 


2  Qi 


Ei{t ) 


where  Qi  =  e0eiflr/ai. 


ICOiQ 


Pi,r(t) 


(6.20) 


(6.21) 


Substituting  the  explicit  forms  of  Ps,q(t )  and  Pt.r(t)  into  the  SVEA  equations,  we  obtain 
the  coupled  SVEA  equations, 


dEs(t) 

dt 

dEiit) 
dt  + 


i(ytq  —  UJso)  + 

%(Qr  ca^o)  V 


jV 

2  Qs 

2Qi 


Es(t) 


Ei{t) 


A^o\/(2) 

4es 

iUioX{2) 

4e* 


exp(zA/cL/2)  sinc(AfcL/2)  E*(t)Ep(t) 


exp(iA/cL/2)  sinc(AfcL/2) 


6.3  Standing  Wave 


To  quantize  a  standing-wave  field, 


Est(z,  t )  =  Est(t)  sin  — — 
we  note  that  the  energy  per  area  stored  in  cavity 

W(t)  —  -  e0eL\Est(t)\2 

where  L  is  the  cavity  length,  is  equal  to  the  energy  per  area  of  the  photon, 


(6.23) 


(6.24) 


W(t)  =  hust\ast(t)\2  (6.25) 

where  u>st  is  the  angular  freuqnecy  of  the  photon  and  ast(t)  is  the  annihilation  operator  of 
the  photon  with  units  of  (photon  number/area)1//2.  The  conversion  from  Est(t )  to  ast(t )  is 
then 


Est®  ~  2V 1  e0esStL  ast (t) 
For  E*t(t )  and  a'st(t),  the  conversion  is 

K(t)  = 


(6.26) 


(6.27) 


6.4  Traveling  Wave 

For  a  traveling  wave  Et(t),  the  power  per  unit  area  is 

m  =  2-|  (6.28) 

zVt 

which  is  equal  to  that  of  the  photon  field, 
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(6.29) 


hit)  =  hut\at(t)\2 

where  at(t )  is  the  annihilation  operator  of  the  photon  held  and  has  units  of  (photon  number/time/area)1/2. 
The  conversion  between  Et(t)  and  at(t)  is  then 

Et(t)  =  \J 2 hx>tr]t  at(t)  (6.30) 

For  E*(t)  and  a\(t),  it  is 

E*(t)  =  \j2hutrit  a\(t)  (6.31) 


6.5  Quantized  coupled  SVEA  (Lamb)  Equations 

The  coupled  SVEA  equations  are 


d  Es(t) 
dt  + 
dEi(t) 
dt 


CJso)  + 
i(Qr  —  ujio)  + 


2  Qs 

2Q, 


Es(t) 


m ) 


UJs  0X(2) 
4es 


exp(iAkL/2 )  sinc(AfcL/2)  E*(t)Ep(t ) 


k^oX(2) 

4e* 


exp(iAkL/2)  sinc(AfcL/2)  E*(t)f£fr{$) 


Using  the  held  normalizations  and  noting  that  rjs ti  =  rj0/nSji,  es>i  =  n2si,  and  c  =  l/e0rj0, 
we  obtain  the  quantized  coupled  SVEA  equations, 

-in  a\(t )  +  VtT 
-in  a\{t )  +  6“(t) 

a«(t)  -  W) 
a/v  <*•(*)  - 
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+ 

rs 

i{£lq  ^so) 

as(t) 

daiit) 

dt 

p; 

i(t2r  —  cujo)  +  — 

Oi(t ) 

rw 

rw 

(6.33) 


where  5“  (f)  and  b("f(t)  are  the  annihilation  operators  of  the  input  and  output  coupling  fields 


(traveling  wave)  at  the  output  mirror,  respectively,  and  have  units  of  (photon  number/area/time)1/2. 


rSji  =  Qq)r/Qs,i  —  As,i[2£,s,i  +  (1  —  7”S)j)]  are  the  total  power  decay  rates  with  ASii  being  the 
free  spectral  ranges,  rS)j  being  the  reflectivities  of  mirrors,  and  being  the  absorption  losses 
of  crystal  per  trip,  and  7 S)j  =  ASij(l  —  rsp)  are  the  power  decay  rates  due  to  the  coupling. 
k  =  kq  exp(fA/cL/2)  sinc(A/cL/2)  with  ac0  =  \  de^c2risrii(tjjsuji)1l2Epi  where  d  =  x^/2 ,  u>s 
and  u>i  are  the  oscillation  frequencies  of  the  signal  and  idler  helds,  respectively,  and  uso  ~  us 
and  Uio  ~  u>i  are  used. 

6.6  Power  Conservation 

We  first  consider  the  case  where  there  is  no  pump,  Ep  =  0,  and  therefore  k  =  0.  The  coupled 
equations  become 


(6.34) 


where  T(  =  i2(Qq  —  ujs0)  +  Ts  and  T'  =  i2(Qr  —  ui0 )  +  T, . 


To  calculate  the  photon  rate  per  area  of  the  signal  held  inside  the  cavity, 


(6.35) 


we  note  that 
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We  hence  have 


d\as(t)  p 
dt 


Is 


—  211  —  T-i  |6i"(()|2-^|6rw|2+fl  -  h  )  (t)+(  1  -  n  )  vf  (i)6“l(i) 


r. 

7s 


(6.37) 


where  we  use  the  commutator  [&*(£'),  feJ(f)]  =  SijS(t'  —  t). 


For  a  lossless  cavity,  Ts  —  j8,  we  then  have 


3M*)I2 

dt 


K(t)\2-\b°ut(t)\2 


(6.38) 


which  says  that  the  photon  rate  per  area  of  the  signal  held  inside  the  cavity  equals  to  the 
difference  of  that  of  the  input  and  output  photons  in  the  same  cavity  mode.  Similarly,  for 
the  idler  held, 


d\oi(t)  |2 
dt 


K(t)\2-\VU\t)\2 


(6.39) 


6.7  With  Pump 


We  next  consider  the  case  where  pump  is  present.  We  note  that 


<*  W- 


dt 


=  al(t) 

=  a\(t ) 


V?(t)  -  y  as(t)  -  in  a\ (t) 


VT,  Kit)  -  h  a,(t) 


(6.40) 


The  photon  rate  per  area  of  the  signal  held  inside  a  lossless  cavity  is  then 
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d\as(t)\- 

dt 


=  Kitfl2  -  l&°ut(f)l'  +  [*«*  as  {t)ai  (t)  -  iK  aj(f)a!(f)] 


(6.41) 


Similarly,  for  the  idler  field,  we  have 


d\ai(t)  p 

dt 


=  ICWI2  -  ir'WI2  +  K  «.(*)<*(*)  -  «  4w4(*)] 


(6.42) 


where  we  use  the  commutator  [a^t'),  Ofc(t)]  =  SjkS(t '  —  t).  Comparing  the  above  two  equa¬ 
tions,  we  see  that 


9k«i-  +  |6jutu)i2  -  icwi2  =  rwr  -  icwr 


<9t 


(6.43) 


6.8  Output  Fields 

The  coupled  SVEA  equations  in  the  time  domain  are 


das[t) 
dt  + 


—i(uso  —  fig)  + 


aa(t ) 


da\(t) 

dt 


i{uiio  flj.)  T 


r* 

2 


<4(f)  +  v'tT  W) 

m  as(t)  +  yfyl  6“f(t) 


(6.44) 


We  will  use  the  Fourier  transform  pairs, 


/(*)  = 


f(5uj)  exp(— iScot)  d&n 


/(M  =  f  /(f)  exp(iSut)  dt  (6.45) 

J  —  OO 

to  obtain  the  coupled  SVEA  equations  in  the  frequency  domain,  where  Su>  is  a  low-frequency 
variable  that  is  relative  to  the  carrier  frequencies  ojs o  or  oji0.  For  the  signal  held,  the  Fourier 
transform  pairs  of  the  annihilation  operator  are 
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oo 


as(t)  =  /  as(u  -  us0)  exp[-i(cu  -  ujs0)t]  d(u;  -  cus0) 

J  —  OO 

1 

as(u  —  uso)  =  as(t)exp[i(u;  -  u;s0)t]  dt  (6.46) 

J  —  OO 

where  u ;  is  the  oscillation  frequency  of  the  signal  held.  For  the  idler  held,  the  Fourier 
transform  pairs  of  the  annihilation  operator  are 


<h(t)  =  /  a^Ui  -  uji0 )  exp[-i((Ui  -  ui0)t]  d(u;»  -  ui0) 

J  —OO 

di(uJi-oJio)  =  —  /  exp[i(cUi  -  uji0)t\  dt 


(6.47) 


where  u>i  =  u>p  —  co  is  the  oscillation  frequency  of  the  idler  held. 


Noting  that 


T 


T  [/*(<)] 

a/m 


-i  POO  C  -i  POO 

—  /  /f(t)exp(i<5u;t)dt  =  |  —  / 


dt 


1 

2tt 


'  —OO 
POO 


9f(t)  ,N1,  1 

—  eXp(i6u*)dt  =  - 


— OO 

°c  f  Q 


f(t)  exp[i(-<5w)t]dt  |  = 


dt 


f(t )  exp (iSuit)  —  )  exp(i&nt)  >  dt 


2vr 


[/(t)  exp(t<5o;t)]^000  —  idujf(Suj)  =  — i5ujf(5uj ) 


(6.48) 


where  we  assume  that  f(t)  is  bounded  so  that  f(t  — >•  ±oo)  =  0,  we  obtain  the  Fourier- 
transformed  coupled  SVEA  equations 


-i(u  -  fig)  as(u>  -  us0)  +  y  as(w  -  u}s0) 

p. 

—  fir)  a\(—uji  +  ujio)  +  —  a\(—uji  +  ojio) 


-in  +  uji0)  +  b™ {u  -  us0) 

in  as( co  -  ujs0 )  +  ^  +  u/J$.49) 
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The  operators  in  the  coupled  equations  are  slowly-varying  envelopes.  To  convert  them  to 
rapidly- varying  analytic  signals,  as(u>),  a* (a;*),  bs(u),  and  6*  (a;*),  we  note  that 


aSti(t)  =  as>i(t)  exp(-iujs0ji0t) 

bs,i(t )  =  bSji(t)  exp(-iujs0ji0t)  (6.50) 


and  therefore 


as,i(u}')  —  +  a;so,io) 

bs,i(u  )  bSji(u)  T  Wso,«o)  (6.51) 

The  coupled  SVEA  equations  for  the  rapidly-varying  analytic  signals  are  then 


r 

-i(uj  -  Qg)  as(u)  +  Y  as(uj)  =  -in  a] (-a;*)  +  y/%  b™(u) 

r  , 

i{ui  -  fir)  a](-uji)  +  y  =  in  as{u)  +  y/%  bf\-Ui) 


(6.52) 


We  rewrite  the  above  equations  in  the  matrix  form, 


IK 


aAu 


-in  i{ui  -  fir)  +  ^  )  \  a\(-uji) 


-  flq)  +  ^ 

—in 

to  obtain  as(co)  and  aj(— ay)  in  terms  of  bf(co)  and  fe‘n^(— ay), 


y/Ts  W 


-OJi 


(6.53) 


f  as(u>) 


(  +  \  1  (  VTs  K(u)  \ 

\  “*«  ~  Mr)  +  %  )  \y/Ti  K\~Ui)  ) 

f  [l{Ui  -  Or)  +  ^\y/Ta  \ 

V _ iKy/Ts _ HO  -  tog)  +  ^y/Ti  )  (  bf(oj) 


~i(UJ  —  fig)  + 


IK 


— in 


f(oy  —  flr)  +  -rf- 
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The  slowly- varying  output  fields,  b°ut(t )  and  b°ut(t),  are  related  to  the  internal  fields,  as(t) 
and  a,i(t),  and  the  input  fields,  bf(t)  and  &‘n(f),  through 


(  <C‘W  \  _(  s/Ts  0  \(  a,(t)  \_(  >C<ti  \ 

l  r*f(«) )  v  °  vit  A  °lw )  v  f >;ntw ) 

Taking  Fourier  transform  of  these  equations  followed  by  convertion  to  rapidly-varying  ana¬ 


lytic  signals,  we  have 


(  V  _(  sFT,  0  \/  a,(w)  W  \  /AM  b(u)  \  / 

Urb-^)/  V  o  fcM 

(6.56) 


=  %)TO  +  5Mr(-^) 

&°utt(-Wi)  =  C'(w)  S“(w)  +  D(u)  S“t(-wi)  (6.57) 

where 

[y  +  *(cai  -  fir)][7s  -  +  i(u  -  fig)]  +  k2 

[^f  -  i(u  -  fi9)]  m  +  i(uJi  -  fir)]  +  k2 

_ _ 

[y  -  i(u  -  fig)][^  +  i(uJi  -  fir)]  -  k2 

_ _ 

[y  -  *(ca  -  fig)][^f  +  i(ui  -  fir)]  -  k2 

[t  ~  ijv  -  flg)][7»  -  T  ~  i(vi  -  fir)]  +  K2  .  , 

[Y  ~  -  toq)][^  +  *(^i  -  fir)]  - 

For  a  lossless  cavity  (ys  =  Ts,  y*  =  Tj),  A(u>),  B(u>),  C(u>),  and  D(u>)  are  related  by  the 
unitary  conditions, 


A(u)  = 
B(u>)  = 
C(u)  = 
D(u>)  = 

and  oji  =  ojp  —  u>. 
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A(u>)C*(u>)  =  B(lo)D*(lo) 


\A(uj)\2-\B{uj)\2  =  1 

\C(u)\2-\D(u)\2  =  1  (6.59) 


6.9  Spatial  Coupling  Factor  and  Gaussian  Beam 

Before  calculating  the  paired  count  rate  and  spectral  Power  spectral  density,  we  will  include 
the  spatial  coupling  factors  at  confocal  condition  into  k. 

The  modified  classical  coupled  SVEA  equations  read 


dEs(t) 
dt  + 
dEiit) 
dt  + 


—  OAjo)  + 
i(Qr  —  Uio)  + 


2  Qs 

2 & 


Es(t) 


Ei(t ) 


~9a 


ius  oX(2) 
4es 


exp(iAkL/2 )  sinc(AA;L/2)  E*(t)Ep(t) 


~ 9% 


k^oX(2) 
4  €i 


exp[iAkL/2 )  sinc(AfcL/2)  E*(t)(Bp§6)) 


where  the  spatial  coupling  factors  are 


2W2W2  2APAZ 

9s  WP2WS2  +  W2W2  +  W2W2  APAS  +  ApAi  +  ASA{ 

2  IP 2  IP 2  2APAS  ,  . 

yi  W2W2  +  W2W2  +  W2W2  APA  +  APAS  +  AsAi  {  ; 

with  IPp,  1PS,  and  IP*  (defined  below)  being  the  1/e2  beam  radii  of  the  pump,  signal,  and 
idler  beams  at  the  focus,  respectively,  and  Ap,  As,  and  Ai  are  the  corresponding  spot  sizes. 
When  the  cofocal  length  equals  to  the  crystal  length,  spot  sizes  are  the  smallest 
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Ap 


Ai 


\pL 

4  np 
A  aL 
4  np 

A iL 

Arii 


(6.62) 


The  field  normalizations  change  accordingly  to 


Ej(t) 


=  2, 


6o€jLAj 


a,j(t)  for  standing  waves 


'  2hUjT]j 

~A~ 


Oj  (t)  for  traveling  waves 


(6.63) 


and  dj(t)  has  unit  of  photon  number  for  standing  waves  and  photon  number/s  for  trvaeling 


waves. 


The  quantized  coupled  SVEA  equations  have  the  same  form, 


das{t) 

dt 

dai(t ) 

dt 


+ 

+ 


i(Clq  ~  Ldso)  +  — 

r,; 

z(f2r  ^io)  T  T 


a>s(t)  =  a\ (t)  +  VtI  W) 

Oi(t)  =  -in  a\(t)  +  ^/7i  &“(*) 


where  k  =  Koexp(iAkL/2)smc(AkL/2)  with 


(6.64) 


7  2  2  (  \  i/2  ip  2\/AsAjAp 

K0  =  2  Ep  •  A2A2  +  A2A2  +  A2A2 


p  l  1  p  S  1  S  I 


(6.65) 


38 


6.10  Paired  Count  Rate  and  Spectral  Power  Density 

We  assume  that  the  total  loss  of  the  signal  and  idler  fields  are  the  same.  The  paired  count 
rate  is  then  equal  to  the  photon  count  rate  of  the  signal  (or  idler)  field  and  given  by 


dwidu2  e<w',e-<W2t  Cb°K*Ki)i>7\u2) 


r*  OO  COO 


’  — oo  J  —  OO 


(6.66) 


Noting  the  commutator  relation,  [6)n(u;i),  ^'(o^)]  =  ^<5^(5 (tui  —  (u2),  we  have 


—  ([^*(wi)  y^(wi)  +  b*(uj  i)  y(wi)][A(w2)  km + b(lj2)  (w2)]) 

=  B*(Wl)B(W2)(i'“(a2l)i'“t(w2)) 

=  B'(uh)B(uh)  T«(wi  -  «2)  (6.67) 

/7T 

The  paired  count  rate  is  then 


Rcav  /  'S'cav  (^ )  d(U 


where  the  spectral  power  density  (spectrum)  of  the  signal  is  given  by 


(6.68) 


U«)  =  T|b(w)|2 

J _ Isli  4  sinc2(A/cL/2) _ ,  . 

2tt  t[4(w  _  p?)2  +  r2][4(tt;.  _  pr)2  +  ff]  -  ik2[4(w  _  p9)(^  _  pr)  +  rsr,f  +  ^ 


In  the  low  gain  regime,  we  drop  the  k2  and  k4  terms  in  denominator.  The  spectrum 
becomes 
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Sca.V  (u  ) 

and  the  paired  count  rate  is 


1  167s ^ li  Ko  sinc2(A/cL/2) 

2vr  [4{u  -  fig)2  +  rf]  [4(07  -  fir)2  +  r2] 


47s7i(f7  +  f-)« o  sinc2(A/cL/2) 

(rs  +  r.^  +  ^-^-a)2 


(6.70) 


(6.71) 


For  rs  ~  Tj  =  T  and  ys  ~  7 «  =  7,  the  spectrum  can  be  written  as 


.  .  1  16  72  Kq  sinc2(A/:;L/2) 

CaV  ^  =  27T  [4(w  -  fig)2  +  f 2]  [4(07  -  flr)2  +  f 2] 

1  16  7 2  Kq  sinc2(AA:L/2) 

27t  {4[(o;  —  ujsq)  +  (a;p  —  fig  —  h2r)/2]2  +  r2}{4[(o;  —  ojsq)  —  (lop  —  fig  —  flr)/2)2]2 

with  the  central  frequency  being  Wso  =  (up  +  fig  —  flr)/2  for  the  signal  held  and  070  = 
Up  —  us0  =  ( ojp  +  flr  —  flg)/2  for  the  idler  held.  The  FWHM  bandwidth  depends  on  the 
detuning  (c op  —  Qq  —  flr)/2 


Awcav  —  \  r2  +  (cap  —  fig  —  fir)2  +  \/2'v/r4  + 


At  exact  phase  matching ,  A k  =  0,  and  zero  detuning ,  5 
paired  photon  rate  has  the  maximum  value  of 


(t Op  —  fig  —  flr)4 

=  (c Op  fig  flr)/2 


(6.73) 
0,  the 


R 


max 


and  the  FWHM  bandwidth  is  minimum 


(rs  + 


(6.74) 


To  hnd  out  the  mirror  rehectivity  r  that  maximizes  the  paired  count  rate,  we  assume  that 
the  signal  and  idler  helds  have  the  same  mirror  rehectivity  rs  =  Ti  =  r  and  crystal  loss 


40 


€s  =  &  =  €■  Then  dRma,Jdr  =  0  gives  r  =  1  -  4f. 


When  phase  matching  (A k  =  0)  but  not  double  resonance  (5  ^  0)  is  satisfied,  the  paired 
photon  rate  can  be  expressed  as 


where  we  assume  Ts  «  Tj  =  T. 


R  = 


R„ 


i  +  i6(h/r)2 


(6.76) 


6.11  Temporal  Second-order  Intensity  Correlation  Function 

The  temporal  second-order  intensity  correlation  function  of  the  signal  and  idler  fields  is  given 

by 


G(2)(ti,t2)  =  (bTt\h)bTt\t{)bT{tl)bT\u 


r*oo  /*oo  /»oo  /»oo 


dcddcnsdccad cc4  eiuJlt2eiuJ2tle~iU3tle~^t2 


'  —  OO  J  —  OO  J  —  OO  J  —  OO 


6”,t(u.i)6r,t(o;2)6r(w3)rt(w4) 


(6.77) 


where  the  expectation  value  can  be  carried  out  by  using  the  commutator  [bm(cOi),  bm(u}j)]  = 
S(u)i  —  LOj)/2TT, 


[C(— o;i)S1sn(— Wi)  +  D(—ui)bi  ^ (coi)]  •  [A* (u2)b1^ (u>2)  +  B*(u2)bll 


•  +  B (a^M]  •  [C*(-^nt(-a;4)  +  D\-^ 


+  C (—oji) B* (u>2) B (uj3)C* (—0J4)  (b^(—uji) (ca2 ) 6)nt (cj3 ) ( — ca4 ) ' 
=  C*(— cJi)24*(a;2)d(a;3)C,*(— ca4)(— )2<5(o;i  +  ta2)(5(a;3  +  cj4) 

Z7T 

+  C(—  (jj\)B*  {u}2)B(u}^)C*  cj4)(— )"h(a;i  —  cc4)h(a;2  —  CJ3) 
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We  then  have 


G<2)(t)  = 


6(t) 


— oo 
2 


A(o;)C,*(o;)e^rda; 

+i?: 


+ 


i 

27T 


|£?(a;)|2du; 


1 

27T 


|C(a;)  |2da; 


?2 

^cav 


(6.79) 


where  r  =  t2  —  ti,  i?2aiJ  is  the  constant  background  level  resulting  from  the  uncorrelated 
photons  in  different  pairs,  and  the  biphoton  wavefunction  4>(t)  in  the  low-gain  regime  (n2 
and  higher-order  terms  are  dropped)  and  a  lossless  cavity  is  given  by 


0(r)  = 


1 

2n 

1 


A{tu)C*(tu)eiu}Tdtu 


— oo 
/•oo 


27T  J — OO  4"  *-^-)(<U  —  Wp  +  —  2"y) 


rdcu 


2a/TJ\  k 


2{up  -nq-  Qr)  +  z(Ts  +  Ti)  \  e*(^-^)T-rir/2  ?  r  >  0 


giVlqT +rsT /2 


,  r  <  0 


(6.80) 


Substituting  0(r)  back  to  G^(t),  we  obtain  the  biphoton  wavepacket  in  the  time  domain 


G«2>(r)-J?L  = 


4r.,r,;  K2  {  eVsT  ,  r  <  0 


(r.  +  r,)2  +  4(Wp  -  n,  -  a r)2  (.  e-r-2  ,  r  >  o 


(6.81) 


6.12  SPDC  with  Quasi-Phase  Matched  Backward- wave 


The  coupled  SVEA  equations  for  the  backward-wave  SPDC  with  quasi-phase  matching  are 


das(co,  z ) 
dz 


in(z)a\(uji,  z)  exp[iAk'(uj)z] 


—iK,(z)as(uj,  z)  exp[— iAk'(uj)z\  (6.82) 

where  Ak'ioj)  =  kp  —  ks(u )  +  k% (ujp  —  u ;)  and  k(z)  switches  its  sign  every  domain,  namely, 
k(z)  =  k o,  —  k0,  Kq,  •••  for  the  first,  second,  third,  ...  domains.  The  domain  length  is  chosen 


da}(uj,z) 

dz 
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to  be  ld  =  Am/2  =  rmr/Km  with  Km  =  kp  -  ks(uj(Q))  +  ki(co^Q)),  and  oj(q)  and  uji(Q)  being 
the  mth-order  quasi-phase  matched  signal  and  idler  frequencies,  respectively. 

In  the  low  gain  regime,  we  can  take  the  driving  fields  on  the  right  to  be  independent  of  z 
and  equal  to  their  values  at  z  =  0  for  the  signal  field  and  z  =  L  for  the  idler  field, 


das(co,  z ) 
dz 

da\(uji,z) 

dz 


iK(z)a\(u}i,  L )  exp[iAk' (cu)z] 

— iK,(z)as(u ,  0)  exp[— iAk'(uj)z] 


(6.83) 


Since  k(z)  is  periodic  in  z  with  period  Am,  we  can  expand  k(z)  as  the  summation  of  its 
Fourier  components, 


where 


k(z  = 


)  =  ^2  cmexp(-i2TTmz/Am) 


(6.84) 


A. 


k(z)  exp(i27rmz / Am)dz 


m  J  o 


1 


A. 


m 

Kq 


/•Am/2 

/■Am 

/  k0  exp(i2irmz/ Am)dz  — 

/  kq  exp(i2irmz/ Am)dz 

Jo 

J  Am/2 

[exp(im7r)  —  1] 


im  7T 

i2tv0/rrm  if  m  is  odd 

0  if  m  is  even 


(6.85) 


By  substituting  the  Fourier  expansion  of  k(z)  back  into  the  coupled  equation,  we  have 


8as(uj,A, )  2H0  ^  1  r  -[a  r  //  \  tv-  I  f  /  rN 

- a - = - >  ,  —  exp{*  Afc  (ca)  -  Km\z}al(ui,  L) 

dz  7T  m 

(m)=—oo 


(6.86) 
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for  the  first  coupled  SVEA  equation,  where  (m)  denotes  that  the  summation  is  only  over 
odd  m.  Direct  integration  from  z  =  0  to  L  then  gives 


as(uj,  z)  =  as(u,  0)  -  —  V] 

7 r  z ' 

(m)=— oo 


exp{i  [Ak'(oj)  —  Km]L} 
mi[Ak'(uj)  —  Km\ 


a\(ui,L) 


(6.87) 


If  Km  is  close  to  A k'(oj),  the  summation  is  dominated  by  the  mth  term.  We  then  have 


cls(uj,z)  =  as(ca,0) - -^-Lexp 

rmr 


iAk(u)L 


sine 


Ak{oj)L 


a\(uJi,L) 


(6.88) 


where  Ak(co)  =  kp  —  ks(co)  +  ki(u>p  —  co)  —  Krn.  Similary,  the  second  coupled  SVEA  equation 


becomes 


aj(o;,0) 


f  /  T  \  0  J 

a'(w,  L) - Lexp 

mir 


-iAk(u)L 


sine 


Ak(u)L 


as{ uJi,  0) 


(6.89) 


By  converting  the  signal  and  idler  held  operators  at  z  =  0  and  L  to  the  rapidly-varying 
analytic  signals  with 


as(u,z )  =  as(co,  z)  exp[iks(uj)z] 

a,i(uji,z)  =  z)  exp[iki(uji)z]  (6.90) 

we  have 

as(cu,L)  =  Ai(u>)as(u>,  0)  +  Bi(u>)a\(uii,  L) 

0)  =  Ci(u;)as(u;,  0)  +  D±(u)a\(ui,  L)  (6.91) 

where 
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A\(u) 

=  exp  [—iks(u)L\ 

B\(u) 

2kq 

= - Lexp 

rmr 

iAk{u)L 

2 

C\{u) 

2^o  , 

= - Lexp 

rmr 

—iAk{u)L 

2 

D\{u) 

=  exp  [-iki{ui)L\ 

sine 


Ak{u)L 


exp{-*[fcs(u;)  +  ki(ui)\L} 


Ak{u)L 

sine 

2 

(6.92) 


The  spectrum  of  the  paired  photons  is  then  given  by 


SiH 


1 

27T 
(Ao)2l2 

\  rnn  / 


2vr 


sine 


Ak{u)L 


where  the  phase  mismatch  can  be  approximated  by 


(6.93) 


Ak(u)  =  kp  —  ks(u)  +  ki{ujp  —  u)  —  Km 

~  kp  /cs(tU(Qj)  T  kiiujp  Kr 

-(u{Q)-Uj) 


n(g)  +  n(8) 


-(W-W{Q)) 


n(g)+  (g) 


(6.94) 


with  rig  and  nf  being  the  group  indices  of  the  signal  and  idler  fields  at  cj(Q)  and  oj^q)  = 
Up  —  U(Q),  respectively.  The  spectrum  can  then  be  written  as 


Si(u)  = 


(Ai i)2  l2 
2v r 


sine 


L(nig)+nJg)) 

2c 


[U  —  U: 


(Q)) 


with  a  FWHM  bandwidth  of 


(6.95) 


Acui  = 


1.77  ire 


(n(g)  n(s)^ 


and  the  paired  count  rate 


(6.96) 
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(6.97) 


Ri 


L 


— oo 


5'i(o;)da; 


6.13  Practical  Considerations 

Six  different  system  designs  are  considered.  The  parameters  are  calculated  based  on  the 
following  equations  and  their  values  are  given  in  Table  1  and  Table  2. 

6.14  Single  Longitudinal-mode 

To  obtain  single  longitudinal-mode,  we  can  only  allow  one  pair  of  signal  and  idler  modes  to 
line  up  within  the  gain  lincwidth.  In  other  words,  the  spacing  between  two  line-up  mode 
pairs  or  the  cluster  spacing  has  to  be  smaller  than  the  gain  linewidth.  The  cluster  spacing 
can  be  estimated  by 


(6.98) 


6.15  Tuning  Rate 


The  signal  and  idler  field’s  cavity  resonance  frequencies  fls  and  fl*  can  be  tuned  by  changing 


the  temperature  T.  The  thermal  expansion  coefficient  of  KTP  is  about  a  =  (dL/dT)/ L  = 
10-5  K_1.  The  tuning  rate  can  be  estimated  as  follows. 


(6.99) 


where  n ^  and  nf  '1  are  the  group  indices  of  the  signal  and  idler  fields. 
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6.16  Mode  Hoping 


As  we  tune  the  signal  and  idler’s  cavity  resonance  frequencies  by  a  total  amount  equal  to 
the  difference  of  their  free  spectral  ranges,  mode-hoping  will  occur.  The  temperature  change 
corresponding  to  a  mode-hop  is 


AT, 


A< 


hop 


d  T 


+ 


df^ 
d  T 


(6.100) 


where  As  and  A.;  are  the  free  spectral  ranges  of  the  signal  and  idler  fields,  respectively. 


6.17  Stabilization 

To  maintain  double  resonance  for  a  mode  pair,  the  total  amount  of  change  of  signal  and 
idler’s  cavity  resonance  frequencies  has  to  be  smaller  than  the  sum  of  signal  and  idler  fields’ 
HWHM  cavity  linewidths.  As  a  result,  the  temperature  has  be  stabilized  to 


where  F  is  the  cavity  finesse. 


A  Tstab 


as  i  A 
2  F  A  2  F 

dOs  |  dOj 
d  T  ~ *■  d  T 


(6.101) 


6.18  Coherence  Length  and  Bandwidth 

The  coherence  length  of  the  biphoton  in  the  time  domain  will  be  equal  to 


Ieoh  = 


1 

F" 


(6.102) 


where  Ts  and  T?;  are  the  power  decay  rates  at  the  signal  and  ideler  frequencies,  respectively. 

If  both  phase  matching  and  double  resonance  are  satisfied,  the  biphoton  bandwidth  could 
be  as  narrow  as 


ACUryiivi 


Fi  +  6UU  +  r?  -  r=  -  r? 


(6.103) 
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6.19  Paired  Count  Rate  and  Spectral  Brightness 

If  both  phase  matching  and  double  resonance  are  satisfied,  for  a  pump  with  power  Ppump  = 
0.01  -^threshold  and  beam  waist  wq,  the  paired  count  rate  could  be  as  high  as 


Rr 


(  \  Z  iZ  4  4  2  2  O 

=  a{ - )  dMe0c 


P 


pump 


mix 


9 


Isli 


7r(wo/V2)2  (rs  +  OCC 


(6.104) 


where  m  is  the  the  order  of  quasi-phase  matching  and  g  ~  0.5.  We  assume  confocal  condition 
and  take  into  account  the  overall  detection  efficiency  a. 


The  spectral  brightness  (paired  count  rate  per  MHz  bandwidth)  will  be  Pm ax/ Acnmin  and 
the  spectral  brightness  normalized  to  mW  pump  power  is  then  Pmax/Ao;minPpump. 

To  the  best  of  my  knowledge,  the  highest  spectral  brightness  (normalized  to  mW  pump 
power)  reported  so  far  for  type-II  phase-matched  biphotons  is  from  Benson  group  in  2009, 
which  is 


Pr, 


Acnmill  Ppump 


330  counts/(s  •  MHz  •  mW)  (894.3  nm  photons) 


(6.105) 


6.20  Without  Cavity 

For  SPDC  without  cavity  enhancement,  the  photon  rate  at  pump  power  Ppump  is 


l292 

I  ns  ±nA 


Pbwd  Cl: 


(6.106) 


where  the  “+”  and  signs  in  the  denominator  correspond  to  backward-wave  and  forward- 
wave,  respectively. 
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Backward 

Wave1 

Cavity-ei 
F  =  1000 

lhanced2 

F  =  1600 

Forward 

Wave3 

F  = 

Pump  Power  =  Threshold  Power/ 100 

8  W5 

766  /iW 

445  fiW 

8  W5 

85 

FWHM  Bandwidth  A/fwhm 

0.081  cm-1 

2.09  MHz 

1.44  MHz 

1.55  cm”1 

2.09 

Biphoton  Length 

~  412  ps 

98  ns 

143  ns 

~  21  ps 

98 

Gain  Linewidth 

0.081  cm-1 

0.081  cm-1 

0.081  cm-1 

1.55  cm”1 

1.55 

Cluster  Spacing 

1.75  cm-1 

1.75  cm”1 

1.75 

Max.  Generation  Rate  (s”1) 

7.10xl06 

1.31xl05 

9.24xl04 

5.42xl09 

3.23x10 

Max.  Detected  Rate  i?max  (s”1)® 

2.09xl04 

385 

272 

1.59xl07 

Rmax/A/FWHM  (s^MHz”1) 

8.6 

184 

189 

342 

Rmax/A/FWHM-Ppump  (s-1MHz-1mW-1) 

0.0011 

240 

424 

0.043 

5328 

Signal  Free  Spectral  Range 

0.094  cm-1 

0.094  cm”1 

0.094 

Idler  Free  Spectral  Range 

0.089  cm-1 

0.089  cm”1 

0.089 

Signal  Resonance  Frequency  Tuning  Rate 

0.15  cm_1/K 

0.15  cm”1/K 

0.15 

Idler  Resonance  Frequency  Tuning  Rate 

0.20  cm_1/K 

0.20  cm”1/K 

0.20 

Mode-hop  Temperature  Change 

13.7  mK 

13.7  mK 

13.7 

Stability  Requirement 

0.3  mK 

0.2  mK 

0.3 

Table  1:  3-cm  crystals. 


Backward 

Wave1 

Cavity-ei 
F  =  1000 

lhanced2 

F  =  1600 

Forward 

Wave3 

Ca 

F  = 

Pump  Power  =  Threshold  Power/ 100 

8  W5 

1  mW 

441  /iW 

8  Ws 

117 

FWHM  Bandwidth  A/Fwhm 

0.24  cm”1 

5.61  MHz 

3.64  MHz 

4.66  cm”1 

5.61 

Biphoton  Length 

~  137  ps 

37  ns 

56  ns 

~  7  ps 

37 

Gain  Linewidth 

0.24  cm”1 

0.24  cm”1 

0.24  cm”1 

4.66  cm”1 

4.66 

Cluster  Spacing 

5.26  cm”1 

5.26  cm”1 

5.26 

Max.  Generation  Rate  (s”1) 

7.10x10® 

2.52x10s 

1.52x10s 

5.42xl09 

6.22x10 

Max.  Detected  Rate  _Rmax  (s”1)6 

2.09xl04 

741 

448 

1.59xl07 

1828 

Rmax/A/FWHM  (s_1MHz  ' ) 

2.9 

132 

123 

114 

326 

Rmax/A/FWHM-Ppump  (s”1MHz”1mW”1) 

0.00036 

126 

279 

0.014 

2795 

Signal  Free  Spectral  Range 

0.28  cm”1 

0.28  cm”1 

0.28 

Idler  Free  Spectral  Range 

0.27  cm”1 

0.27  cm”1 

0.27 

Signal  Resonance  Frequency  Tuning  Rate 

0.15  cm”1/K 

0.15  cm”1/K 

0.15  cm 

Idler  Resonance  Frequency  Tuning  Rate 

0.20  cm”1/K 

0.20  cm”1/K 

0.20  cm 

Mode-hop  Temperature  Change 

41.1  mK 

41.1  mK 

41.1 

Stability  Requirement 

0.8  mK 

0.5  mK 

0.8 

Table  2:  1-cm  crystals. 
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